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ABSTRACT 

We assess the impact of Galactic synchrotron foreground removal on the observa- 
tion of high-redshift quasar HII regions in redshifted 21-cm emission. We consider the 
case where a quasar is observed in an intergalactic medium (IGM) whose ionisation 
structure evolves slowly relative to the light crossing time of the HII region, as well as 
the case where the evolution is rapid. The latter case is expected towards the end of the 
reionisation era where the highest redshift luminous quasars will be observed. In the 
absence of foregrounds the fraction of neutral hydrogen in the IGM could be measured 
directly from the contrast between the HII region and surrounding IGM. However, we 
find that foreground removal lowers the observed contrast between the HII region and 
the IGM. This indicates that measurement of the neutral fraction would require mod- 
elling to correct for this systematic effect. On the other hand, foreground removal does 
not modify the most prominent features of the 21-cm maps. Using a simple algorithm 
we demonstrate that measurements of the size and shape of observed HII regions will 
not be affected by continuum foreground removal. Moreover, measurements of these 
quantities will not be adversely affected by the presence of a rapidly evolving IGM. 

Key words: cosmology: diffuse radiation, large-scale structure, theory - galaxies: 
high redshift, intergalactic medium 



1 INTRODUCTION 

The reionisation of cosmic hydrogen by the first stars and 
galaxies was an important milestone in the history of the 
Universe (e.g. |Barkana fc Loeb 20011. A powerful tool for 



(e.g. |Tozzi et al.||2000| |Furlanetto et al.||2004| |Loeb fc Zal 



study of the reionisation history will be provided by the 
redshifted 21-cm emission from neutral hydrogen in the in- 
tergalactic medium (IGM), and several probes of the reioni- 
sation era in redshifted 21-cm emission have been suggested. 
These include observation of the emission as a function of 
redshift averaged over a large area of sky. This observation 
would provide a direct probe of the evolution in the neu- 
tral fraction of the IGM, and is referred to as the global 



step ( |Shaver et al.|1999| |Gnedin fc Shaver|2004| |Furlanetto 



20061. Unless reionisation is uniform throughout the whole 



darriaga|2004||Wyithe fc Morales|2007||liiev et al.|2006| >. Fi- 

nally, observation of individual HII regions will probe quasar 
physics as well as the evolution of the neutral gas ( |Wyithe 
fc Loeb|2004||Kohler et al.|2005|[Valdes et al.|2006| >. |Kohier 
et al.| ( |2005[ ) have generated synthetic spectra using cosmo- 
logical simulations and conclude that quasar HII regions will 
provide the most prominent individual cosmological signals. 
Work by |Datta et al.| ( |2007[ > has focused on the detection 
of ionised bubbles in redshifted 21-cm maps. Their results 
suggest it may be possible to blindly detect spherical HII 
regions of radius J> 22 comoving Mpc during the epoch of 
reionisation. More recently, Geil & Wyithe ( 2007 1 (hereafter 
GW07 1 have studied the impact of a percolating IGM on the 



IGM however, the global step will be very difficult to detect 
beneath the bright foreground. A more powerful probe will 
be provided by observation of the 21-cm power spectrum of 
fluctuations together with its evolution with redshift. This 
observation would trace the evolution of neutral gas with 
redshift as well as the topology of the reionisation process 



detection of HII regions, and shown that quasars will leave 
a detectable imprint until very late in the reionisation era. 

Various experiments are planned to measure 21-cm 
emission from the pre-reionisation IGM, including the Low 
Frequency ArrajQ (LOFAR) and the Murchison Widefield 
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Arra>[^] (MWA). In addition to their physical configuration, 
the sensitivity of these telescopes to the epoch of reionisation 
will be limited by the difficulty of achieving a perfect calibra- 
tion and by problems associated with foregrounds and the 
ionosphere. There has been significant discussion in the lit- 
erature regarding low-frequency foregrounds. While fluctua- 
tions due to foregrounds are expected to be orders of magni- 



tude larger than the reionisation signal (e.g. Di Matteo et al. 
2002 Oh & Mack 2003 ) , foreground spectra are anticipated 
to be smooth. Since the reionisation signature includes fluc- 
tuations in frequency as well as angle, it has therefore been 
proposed that they be removed through continuum subtrac- 
tion (e.g. |Gnedin fc Shaver||2004| |Wang et al.||2006[ ), or us- 
ing the differences in symmetry from power spectra analysis 
( |Morales fc Hewitt|2004}|Zaldarriaga et al.|2004| . 

In this paper we focus on the observation of individual 
HII regions around known quasars, with emphasis on the ca- 
pabilities of the MWA. In particular we address the removal 
of a bright Galactic foreground from the redshifted 21-cm 
signal. Our analysis extends the model and ideas presented 
in |GW07| by including evolution of the percolating IGM. 
We begin by summarising our semi-numerical model for the 
ionisation structure of the IGM (§[5]). We show examples of 
model HII regions in §[3] (both with and without evolution 
of the IGM). We then discuss the influence of foregrounds 
in §[4] and describe the observation of the HII region with a 
low-frequency array in §[5] including discussion of foreground 
removal and reconstruction of the HII region shape. The ef- 
fects of different foreground removal models on HII region 
observables are investigated in §[6] We present our conclu- 
sions in §[7] Throughout we adopt the set of cosmological 
parameters determined by WMAP (Spergel et al. 20071 for 
a flat ACDM universe. 



2 SEMI-NUMERICAL MODEL FOR THE 

GROWTH OF QUASAR HII REGIONS IN A 
PERCOLATING IGM 

In this section we describe our model for the reionisation 
of a three-dimensional volume of the IGM using a semi- 
numerical scheme (|Geil fc Wyithe"|p007l |Zahn et al.]|2007} 
Mesinger fc Furlanetto|2007[ ). We follow the implementation 
described in |GW07[ and refer the reader to that paper for 
details of the model. The model employs a semi-analytic 
prescription for the reionisation process, which is combined 
with a realisation of the density field. Each of these aspects 
is summarised briefly below. 



2.1 Density-dependent model of global 
reionisation 

The semi-analytic component of our model computes the 
relation between the local dark matter overdensity and the 
reionisation of the IGM, and is based on the model described 



by |Wyithe fc Loeb| ( |2007| ) and |Wyithe fc Morales| (20071 
The model includes the overabundance of galaxies in over 



dense regions that results from galaxy bias ( |Mo fc White 
1996), as well as the increase in the recombination rate in 



overdense regions, and ionisation feedback which suppresses 
low-mass galaxy formation in ionised regions. This model 
predicts the sum of astrophysical effects to be dominated 
by galaxy bias, and that as a result, overdense regions are 
reionised first. This leads to the growth of HII regions via 
a phase of percolation during which individual HII regions 
overlap around clustered sources in overdense regions of the 
universe. The output of the semi-analytic model is evolution 
of the ionisation fraction by mass Qs.r of a particular region 
of scale R with overdensity 8. 

Throughout this paper we consider a model in which the 
mean IGM is reionised at z — 6 (Fan et al. 112006 Gnedin & 



|Fan|20 06 Whit e et al.|2003| |. We assume that star formation 
proceeds in halos above the hydrogen cooling threshold in 
neutral regions of IGM. In ionised regions of the IGM star 
formation is assumed to be suppressed by radiative feedback. 



2.2 The ionisation field 

We construct an ionisation field based on a Gaussian random 
field for the matter overdensity, combined with the value of 
the ionised fraction Qs,r as a function of overdensity and 
smoothing scale. We repeatedly filter the linear density field 
at logarithmic intervals on scales comparable to the box size 
down to the grid scale size. For all filter scales, the ionisa- 
tion state of each grid position is determined using Qs,r and 
deemed to be fully ionised if Qs.r ^ 1. All voxels within a 
sphere of radius R centered on these positions are flagged 
and assigned Qs,r ~ 1, while the remaining non-ionised vox- 
els are assigned an ionised fraction of Qs,r { min , where i?f im i n 
corresponds to the smallest smoothing scale. A voxel forms 
part of an HII region if Qs,r > 1 on any scale R. In this paper 
we present simulations corresponding to a linear density field 
of resolution 256 3 , with a comoving side length of 512 Mpc. 
We label positions within the simulation with the vector 
x = (xx, X2, X3), and define the center of the box as the ori- 
gin. We take the a;3-axis to be parallel to the line of sight, 
and show results using dimensions of proper distance. The 
line-of-sight coordinate (X3) can be equivalently described 
using redshift or frequency. We have chosen the side length 
of 512 comoving Mpc to match the 32 MHz bandpass of the 
MWA, centered on a quasar at 2 = 6.65. 



2.3 The evolution of quasar HII regions in a 
percolating IGM 

A semi-numerical simulation generates a realisation of the 
three-dimensional ionisation field at a particular value of 



proper time (GW07). However, late in the reionisation era, 



http: / / www.haystack.mit.edu/ast / arrays/mwa/ 



where we expect to observe the high-redshift quasars, the 
light travel time across a quasar HII region may be compara- 
ble to the timescale over which the IGM evolves significantly. 
This evolution may have a negative impact on the observ- 
ability of quasar HII regions. In this paper we therefore de- 
scribe semi-numerical calculations of an ionisation field for 
two different cases of IGM evolution. In the first we assume 
a non-evolving IGM with a slowly evolving quasar. In the 
second we consider a model where the line-of-sight axis in- 
cludes the finite light travel time, and evolution of the IGM 
ionisation and quasar luminosity. The first case corresponds 
to the model in lGW07 l The additions to this model that are 
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required for the second case (with an evolving IGM) are as 
follows. 

Unlike an iV-body simulation which evolves an ionisa- 
tion field in time, a semi- numerical calculation computes the 
ionisation field at a single instant in time. To include the 
finite speed of light we must therefore compute many sim- 
ulations at closely separated redshifts and stack slices with 
line-of-sight distances such that the photons from all slices 
reach the observer at the same time. The redshift of a slice 
is 



ZQ + 



celt 
dz 



(x 3 - xl), 



(1) 



where X3 is the proper distance coordinate along the line of 
sight, and X® is proper the distance to redshift zq . 



2.4 Inclusion of quasars in the semi-numerical 
scheme 

In this section we describe the effect of the quasar on the 
evolving ionisation structure of the IGM. The quasar has 
an influence which evolves with time. We assume the quasar 
turns on(off) at a redshift z on ( ff) (corresponding to a proper 
time ion(off))) an d has a lifetime t\ t — t Q g — t on . Following 
the calculation in |GW07| we parameterise the influence of 
the quasar using the distance R q , which is the radius of the 
region that would have been ionised by the quasar alone in 
a fully neutral IGM 







-Rq — 



-Rq.tot (t/ton) 



1/3 



R, 



-q,tot 



Z ^> Z Q n 

ZoS < Z < Z 

Z < Z S, 



(2) 



where i? q ,tot is the maximum radius of the HII region that 
would have been generated by the quasar in a fully neutral 
IGM. 

To include the influence of a quasar in an evolving IGM 
we generate an ionisation cube for the IGM as a function of 
look-back time using the procedure described above. How- 
ever in doing so we use equation (J3| to include the time- 
dependent quasar contribution to the ionising photon bud- 
get. In a region of radius R containing a quasar, the cumu- 
lative number of ionisations per baryon will be larger than 
predicted by the semi-analytic model. We include quasars 
in our scheme by first computing the fraction of the IGM 
within a region of radius R centered on the position x that 
has already been ionised by stars. We then add an additional 
ionisation fraction equal to the quasar's contribution 



i(x) = 



Rq 



(3) 



where x q is the position of the quasar, and R q is computed 
via equation |2j|. In contrast to stellar ionisation the quasar 
contribution comes from a point source, and so the contri- 
bution to Q originates from a single voxel only (rather than 
all voxels within |x — x q |). Following this addition we filter 
the ionisation field as described in § |2.2| 

The value of i? q ,tot is subject to large uncertainties, 
including the quasar duty cycle, lifetime and luminosity. 
Assuming a line-of-sight ionising photon emission rate of 
N y /4n photons per second per steradian, a uniform IGM 
and neglecting recombinations, the comoving line-of-sight 



extent of a quasar's HII region (from the quasar to the front 
of the region) observed at time t is ( White et al.|2003 1 

1/3 



R, 



q.tot 



30Mpc x~l /3 



10 57 10 7 yr 



1 + z 
7.5 



(4) 



where xm is the neutral fraction. In this paper we show 
examples of i? q ,tot ~ 34 comoving Mpc, consistent with the 
lower limit around the luminous SDSS quasars at z > 6 (Fan 
et al.|2006 l. 



3 EXAMPLES OF QUASAR HII REGIONS 

We compute the three-dimensional ionisation structure for 
our reionisation model using the prescription outlined in §[2] 
The contrast in 21-cm brightness temperature corresponding 
to a region in the simulation with ionisation fraction Q and 
overdensity S is 



T(S,R) = 22 mK (1-Qs,r) 



7.5 



1/2 



(1 + 5). (5) 



Here we have assumed that the spin-temperature of neutral 
hydrogen in the IGM is much larger than the cosmic mi- 



crowave background (CMB) temperature (Ciardi & Madau 
|2003[ ), and ignored the enhancement of the brightness tem- 
perature fluctuations due to peculiar velocities in overdense 
regions ( |Bharadwaj fc Ali|2005| |Barkana fc Loeb|2005| ). Be- 
cause the ionisation state of the IGM around the highest 
redshift quasars is not known, we consider the two cases de- 
scribed in § |2.3| which bracket the ionisation conditions sur- 
rounding a quasar prior to the completion of reionisation. 

Resulting 21-cm brightness temperature maps are pre- 
sented in Figures [TJ]3] In the upper row of Figure [I] we show 
snapshots of the model brightness temperature T mo d c i for 
slices through a realisation at a time when the mean neu- 
tral fraction of the IGM is w 15%, corresponding to redshift 
z — 6.65 (assuming Rq = iZq.tot). This simulation was com- 
puted at fixed proper time, and so represents the model with 
slow evolution in the reionisation process. Each slice is 3 co- 
moving Mpc deep. We show three aspects of the HII region, 
with slices viewed from the front, side and top. These slices 
are similar because the simulation is statistically isotropic. 
The position of the quasar within the simulation cube is 
plotted as a square point. 

In Figures [2] and [3] we show results for the model in 
which finite light-travel time is included, showing the effect 
of rapid evolution in the ionisation state of the IGM towards 
the end of the reionisation era. In this simulation the center 
of the box is at z = 6.65, at which time the IGM has a 
neutral fraction of w 15%. The quasar is assumed to have 
■Rq.tot = 34 comoving Mpc, a lifetime of 2 x 10 7 yr and be 
centered on z — 6.65. The position of the quasar within 
the simulation cube is plotted as a square point. The upper 
row of Figure [2] shows three aspects of the HII region in 
this model, with slices viewed from the front, side and top. 
The evolution of the IGM is clearly seen in this figure, with 
the percolation process completing between the "back" of 
the box [xa > 0) and the "front" of the box (2:3 < 0). In 
Figure [3] we show a sequence of slices between z — 6.74 and 
6.56. These slices show the thinning of neutral gas towards 
lower redshifts. The presence of an HII region around the 
quasar can be seen towards the middle of the redshift range. 
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Front view Side view Top view 




-20 20 -20 20 -20 20 

Xl (Mpc) x 3 (Mpc) Xi (Mpc) 



Figure 1. A quasar HII region in a non-evolving IGM. Three aspects of the HII region are shown, with slices through the center of the 
box when viewed from the front, top and side. Each slice is 6 Mpc thick, which corresponds to ~ 3 MHz along the a;3-axis. In observed 
units, the cube is ~ 3.3 degrees on a side and 33 MHz deep. The model, foreground plus model, and observed maps following foreground 
removal are shown in the upper, central and lower panels respectively. The shape of the HII region extracted from the foreground removed 
cubes using the method described in the text are also plotted. The quasar was assumed to have R q ,tot = 34 comoving Mpc, and its 
position is plotted as the square point. The mass-averaged IGM neutral fraction is assumed to be 15%. For the purposes of display we 
have used the central redshift of the box to convert between proper and comoving distances throughout this paper. 



4 THE CONTINUUM FOREGROUND 

Foreground contamination and its removal could have sig- 
nificant consequences for the detectability of the weak, red- 
shifted 21-cm signal. Indeed foreground contamination will 
be brighter than the cosmological 21-cm signal by 4-5 or- 
ders of magnitude. The three main sources of foreground 
contamination of the 21-cm signal are Galactic synchrotron 
(which comprises ~ 70%), extragalactic point sources (27%) 
and Galactic bremsstrahlung (1%) ( Shaver et al.|1999 |. The 
frequency dependence of these foregrounds can be approxi- 



mated by power laws with a running spectral index ( Shaver 
et al.||1999| |Tegmark et al.||2000[ ). While the sum of power 
laws is not in general a power law, over a relatively nar- 
row frequency range (such as that considered in this paper 
where Av/v < 1), a Taylor expansion around a power law 
can be used to describe the spectral shape. We therefore also 
approximate the sum of foregrounds as a power law with a 



running spectral index, and specialise to the case of Galac- 
tic synchrotron emission, which dominates the foreground^] 
For more detailed foreground models, including synchrotron 
emission from discrete sources such as supernova remnants, 
and free- free emission from diffuse ionised gas, sec Jclic et al. 



3 The primary effect of other foregrounds will be to slightly alter 
the spectral index and the frequency dependence of the spectral 
index, and also (in the case of unresolved radio point sources) the 
angular structure of the foreground at small scales. Neither of 
these effects is important as our foreground removal technique is 
not sensitive to the specific value of the spectral index we adopt. 
As for angular fluctuations, these correspond to zero-point fluc- 
tuations along different lines of sight, which are immediately re- 
moved by foreground cleaning. We have experimented with in- 
creasing angular fluctuations by a factor of 10, and found little 
difference. 
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Front view Side view 




[Mpc) X3 (Mpc) (Mpc) 

Figure 2. As per Figure [l] this time for a quasar in an evolving IGM. The quasar was assumed to have R q ,tot = 34 comoving Mpc 
and a lifetime of 4 X 10 7 yr centered on z = 6.65, which is also the redshift at the center of the simulation box. The mass-averaged IGM 
neutral fraction at the redshift of the quasar is ?s 15%. 



(20081. We assume that the brightest point sources have 



been removed and only consider unpolarised foregrounds. 

The intensity of Galactic synchrotron emission varies as 
a function of both sky position and frequency. We model the 
frequency and angular dependence of Galactic synchrotron 
foreground emission as follows. We first construct a realisa- 
tion of the angular fluctuations in the foreground (around 
the mean brightness temperature) at a particular frequency 
fo using the relation 



rcijvo) 

2tt 



2-/3 



(6) 



where the frequency dependence of the fluctuation ampli- 
tude is given by 



msyn 

1 lo 



A syn 
<0 



150 MHz 



and the constants have values a s 

and A= yn 
Ui 



Wang et al. 2006 1. The latter two values are extrapolatec 



-Aa syn logioH/(150MHz)] 



(7) 

= 2.55, Aa syn = 0.1 



25 K ( [Shaver et al.|i"999} [Tegmark et al.| [ 



2000 



from 30 GHz CMB observations ( McQuinn et al.|2006 1 since 



at present there are no data at lower frequencies with arc 
minute resolutiorj^] We then add a mean sky brightness 
T syn = 165(f /185MHz)- 2 ' 6 K to our two-dimensional re- 
alisation of brightness temperature fluctuations AT syn (8) 
(although interferometers will generally not be sensitive to 
the temperature zero-point). This foreground plane is then 
extended into three dimensions by extrapolation, using the 



4 Although the high frequencies of CMB experiments make 
them less than ideal for understanding diffuse Galactic emis- 
sion at frequencies relevent to high-redshift 21-cm tomogra- 
phy, several publicly available unpolarised, large area radio sur- 
veys have been compiled to form a Global Sky Model (see e.g. 
http://space.mit.edu/~angelica/gsm/index.html). This, together 
with the 843 MHz sub-degree resolution and ~ ljybeam -1 rms 
sensitivity SUMMS survey ( |Bock et al.|1999| for 6 < -30°, will 
provide a better practical continuum foreground model for the 
MWA target fields (Randall Wayth, private communication). 
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z=6,69 



-20 

K 1 (Mpe) 



z=6.71 
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iti (Mpc) 



x 1 (Mpc) 



Figure 3. Slices through the simulation of a quasar in an evolving IGM for 6.56 z $J 6.74, viewed from the front. The shapes of the 
HII regions extracted from foreground removed cubes using the method described in the text are also plotted. 



running power law form along each line of sight 



\e,u) 



T syn + AT s; 



-Aa ayn logi [i//i/ ] 



(8) 



In this paper we have not modeled the spatial variation of 
a and Aa syn , although the effect of this is negligible, as 
we have verified in other work (Petrovic & Oh, 2008, in 
preparation). Most of the error in foreground fitting comes 
from "noise" due to the instrument and the cosmological 
21-cm signal. 

The total observed brightness temperature becomes 
Tobs = T syn + T mo dci- Maps of the continuum foreground 
model plus 21-cm signal are shown in the second row of pan- 
els in Figures [l] and [2] The foreground completely swamps 
the HII region signal, which is not visible in the combined 
maps. The mean and variance of the foreground brightness 
(at the central frequency of these cubes) are 1.7 x 10 5 mK 
and 8.7 x 10 3 mK respectively, which should be compared 
to the ~ 10 mK 21-cm fluctuations. The front, side and top 
views are shown as before. The front view shows the an- 
gular structure of the foreground, while the side and top 
views show the coherent structure along the line of sight, 
owing to the smooth frequency dependence of the assumed 
foreground emission. 



5 SYNTHETIC OBSERVATIONS OF MODEL 
HII REGIONS 

In this section we discuss the response of a low-frequency 
array to the model HII regions described in §|3] including 
foreground subtraction. We then estimate the potential for 
measurement of HII region properties, including the HII re- 
gion shape and the global neutral fraction. 

5.1 Fundamental sensitivity limit to the 21-cm 
signal 

To compute the response of a low-frequency array to the 
model HII region we follow the prescription given in |GW07] 
Radio interferometers measure a frequency-dependent, com- 
plex visibility V(\J,v) (Jy) for each frequency channel and 
baseline U in their configuration. The measured visibility 
is, in general, a linear combination of signal, foreground and 
system noise, 



V(U, u) = V S (U, v) + V F (U, v) + V N (U, u), 



(9) 



where Vs is the signal, Vf the contribution due to foreground 
sources and Vn the system noise. The root mean square of 
noise fluctuation per visibility per frequency channel is given 
by 

2k-BT sys 



AVn = 



Aesy/tv Al 



(10) 
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where T sys is the system temperature (K), A e g the effective 
area of one antenna (m 2 ), t\j is the integration time for that 
visibility (s), Av is the frequency bin width (Hz) and fee is 
the Boltzmann constant. Instrumental noise is uncorrelated 
in the frequency domain and is assumed to be Gaussian. 
As shown by McQuinn et al. (2006), the average integration 



time t\j that an array observes the visibility U is 



tv 



■n(U), 



(11) 



where ti nt is the total integration time and n(U) is the 
number density of baselines that can observe the visibil- 
ity U. The transverse wavenumber k± is given by k± — 
27r|U|/r cm (2;), where r cm (z) is the proper distance to the 
point of emission. 

We present simulated results for the MWA, which will 
consist of 512 tiles, each with 16 cross-dipoles with 1.07m 
spacing and an effective area A e ^ « 16(A 2 /4) for A <J 2.1m. 
The system temperature at v < 200 MHz is sky dominated 
and has a value T sys ~ 250 [(1 + z)/7] 2 ' 6 K. This value is an 
estimate of the average sky brightness due to Galactic syn- 
chrotron emission over the primary beam. Following Bow- 
man et al. ( 2006 1, we assume a smooth antenna density pro- 
file pant oc r~ within a 750 m radius. For the MWA, the 
aperture is filled within a radius of 18 m. In order to main- 
tain high sensitivity we truncate the naturally weighted vis- 
ibility data using a filled circular aperture of radius a with 
a corresponding beamwidth (full width at half maximum) 
6h = 0.705l7mi x , where (7 m ax = a/A. In this paper we show 
results computed for # D = 3.2', and assume an integration 
time of tint = 100 hours. 

We simulate the thermal noise in a three-dimensional 
visibility- frequency cube using equations (10 1 and (111, al- 



lowing us to account for a particular array configuration and 
observation strategy. We then perform a two-dimensional in- 
verse Fourier transform in the uw-plane for each frequency 
channel in the bandwidth, which gives a realisation of the 
system noise in the image cube (i.e. sky coordinates). Us- 
ing the linearity of both the Fourier transform and equation 
JjjJ, the observed specific intensity can then be found ei- 
ther by first adding visibility terms and then transforming 
to image space, or by adding the realisation of each visibility 
component in image space. Finally, we apply the Rayleigh- 
Jeans approximation to express the observed flux in terms of 



brightness temperature T using dB v /dT 
B v is the specific intensity. 



2/cb/A , where 



5.2 Continuum foreground subtraction 

It has been suggested that spectrally smooth foregrounds 



could be removed through continuum subtraction ([Gnedin 
& Shaver|2004||Wang et al.|2006||Jelic et al.|2008h. We model 



Table 1. Parameters describing the HII region shape in models 
of non-evolving and evolving IGM, corresponding to Figures [l]and 
[2] respectively. The values in parentheses are the offset relative to 
the volume-averaged radius (i? 3 ) 1 / 3 . 





evolving 


; IGM 


non-evolving IGM 




FGR map 


model 


FGR map 


model 




7.1 (-1.7) 


7.0 (-1.6) 


7.1 (-0.2) 


6.8 (-0.6) 


A 2 


9.0 (+0.2) 


8.4 (-0.2) 


7.6 (+0.1) 


7.6 (+0.2) 


Bi 


7.1 (-1.7) 


7.4 (-1.2) 


7.3 (-0.1) 


7.3 (-0.1) 


B 2 


10.3 (+1.5) 


10.4 (+1.8) 


7.6 (+0.2) 


7.5 (+0.1) 


c\ 


11.0 (+2.2) 


10.5 (+1.9) 


7.1 (-0.3) 


7.4 (0.0) 


c 2 


5.6 (-3.2) 


5.6 (-3.0) 


7.6 (+0.5) 


7.6 (+0.2) 


(R 3 ) 1 / 3 


8.8 


8.6 


7.4 


7.4 



where v is the observed frequencjj^] We begin with results 
for ni og = 3 (which we label model II) , but return to discuss 
different foreground removal models in § [6] 

We perform a fit of this functional form along the 
line of sight to each spatial pixel in the simulation cube. 
We then subtract the best-fit, leaving residual fluctuations 
around the foreground emission. The resulting foreground 
subtracted maps are plotted in the lower rows of Figures [l] 
and [2] Inspection of these maps indicates that the main 
features of the HII regions are still present. However, fore- 
ground removal has significantly lowered the contrast of the 
maps relative to the input model. This is because foreground 
removal erases modes in the signal with wavelengths com- 
parable to the bandpass, leaving a foreground-cleaned data 
cube composed of residuals, with a near zero- mean bright- 
ness temperature. 



5.3 The shape of HII regions 

In our model we have assumed isotropic quasar emission, 
which is probably not correct. However even in this case 
the resulting HII region may not be spherical owing to the 
inhomogeneous nature of the stellar reionisation around the 
quasar HII region, as well as the effects of evolution which 
could serve to elongate the HII region along the line of sight 
HWyithe et al.|2005||Yu|2005) . 

In an inhomogeneously ionised IGM the meaning of 
"shape" is not clear. For the purpose of illustration, we de- 
fine the shape of the HII region to be the shape of a spheroid 
that maximises the contrast in 21-cm emission, 



-^contrast — (Tout) - (Tin), 



(13) 



between two concentric shells, where (Tout) and (Ti n ) are the 
mean emission within shells of thickness AR outside and in- 
side a spheroidal surface respectively. We then parameterise 
this spheroid by 



+ 



+ 



= 1, 



(14) 



the continuum using a polynomial of the form 



where a = 01 for xi < and a = a 2 for x\ > 0, b = 61 



lo: 



gio(r fg ) = ^a[io glo Mr 



5 Note that when compared with the frequency evolution as- 
(12) sumed for our synchrotron model, the coefficients C\ and C 2 

correspond to a sy n and Aa 3yn respectively. 
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Table 2. Values of recovered neutral fraction for the four fore- 
ground models. The model values are listed in parentheses. 



model 


index 


evolving IGM 


non-evolving IGM 


I 


"log = 2 


0.07 (0.15) 


0.10 (0.14) 


11 


"log = 3 


0.07 (0.15) 


0.10 (0.14) 


ill 


"log = 4 


0.04 (0.15) 


0.04 (0.14) 


IV 


"lin = 4 


0.04 (0.15) 


0.05 (0.14) 



for X2 < and b = b 2 for x 2 > 0, and c = Ci for x 3 < 
and c = C2 for X3 > 0. To estimate the shape of the HII 
region we therefore maximise the quantity T con trast by vary- 
ing the parameter set p = (ai, <22, &i, b%, Ci, C2). We assume 
AR = 1 Mpc and that the spheroid is centered on the quasar 
(we also assume the position and redshift of the quasar has 
been accurately measured by optical/IR observations). In 
order to evaluate the prospect of measuring the shape of 
a quasar HII region, this maximisation is performed in the 
noisy, foreground subtracted data sets shown in Figures [l] 
and[l 

The resulting optimised spheroids are plotted over both 
model and foreground subtracted noisy maps in Figures [I] 
and [2] The numerical values of the fitted parameters p are 
listed in Table [l] Although the maximisations were per- 
formed in the foreground removed maps, the shapes of the 
recovered HII regions faithfully represent the input model. 
For example, in the evolving IGM model (Figure [5| the fit- 
ted spheroid describes the extension of the observed HII re- 
gion towards the observer along the (negative) a^-direction. 
This faithful reproduction holds for both the non-evolving 
and evolving models, and indicates that foreground removal 
does not erase information regarding HII region shape. To 
quantify this statement we also perform the maximisation 
of Tcontrast on the corresponding input models and list the 
resulting parameters p in Table[l]for comparison. The shape 
of the HII region, including the asymmetries, are very simi- 
lar in the two cases. We also calculate the volume-averaged 
radii (i? 3 ) 1//3 and list these values in Table [l] The values in 
parentheses are the offset of the parameters p with respect 
to this volume-averaged radius. These offsets are similar be- 
tween the input model and noisy foreground maps, reaffirm- 
ing that the foreground removal process does not degrade 
the ability to measure HII region shape and size. 



5.4 Estimate of neutral fraction 

Based on the fact that the HII region shape is not known 
a priori, [GW07| investigated the constraints on neutral frac- 
tion that could be placed using a narrow beam along the 
quasar line of sight. However, in the previous section we 
have argued that the spheroidal shape of the HII region can 
be accurately reproduced from foreground subtracted maps. 
We can therefore use the full data cube to determine the con- 
trast of emission that originates from regions that are inside 
and outside of the HII region respectively, rather than just 
along one line of sight. To achieve this we fit planes of 21-cm 
emission 



T in (x) = Tin + A in X! + B m X 2 + Ci„ X 3 



(15) 



and 

Tout(x) = Tout + Amt X! + Bout X 2 + Cout X S (16) 

to regions of the simulations of foreground subtracted 21-cm 
intensity T(x) = T b s (x) — Tf g (x) that are inside and outside 
the HII region respectively. This method was adopted in 
order to fit for both the mean value and evolution of global 
neutral fraction. Using these solutions we can estimate the 
global mass-averaged neutral fraction at the redshift of the 
quasar, which is obtained from 

Tout - T in fl + z\- 1/2 

Xm = ^2rnR- v TTJ (17) 

Performing these fits for the models shown in Figur es [T] and 
[5] gives xm = 0.10 and 0.07 respectively (see Table |2|. Ide- 
ally, the coefficients A in / out and -B ln / ou t should vanish. The 
value of AC = C ou t — Cout could in principle be used to mea- 
sure the evolution in global neutral fraction, however noise 
and foreground removal will prevent this measurement in 
practice. 

For comparison, we have also computed the correspond- 
ing neutral fractions from fits to the input model, which 
yield values of :ehi = 0.14 and 0.15 respectively. Thus, fore- 
ground removal erases a significant fraction of the contrast 
between the HII region and the IGM, which quantifies the 
observations made from inspection of Figures [l]and[2] m ear- 
lier sections. We note that the fraction of contrast that is 
lost differs between these two examples, owing to the large 
difference in evolutionary properties of the IGM. This sug- 
gests that the effect of foreground removal on the inferred 
neutral fraction could only be reliably estimated via detailed 
numerical modelling. 



6 SENSITIVITY TO THE FOREGROUND 
REMOVAL MODEL 

Before concluding, we investigate the effect of different 
foreground removal models on the HII region observables. 
Specifically, in addition to model II (ni og = 3), we choose 
polynomials in \og 10 (v) with ni og = 2 (labelled model I, 
corresponding to the form of the assumed synchrotron fore- 
ground) and nio g = 4 (labelled model III). We also consider 
a polynomial in v [as opposed to log 10 (^)] of degree nr, n = 4 
(labelled model IV). 

We find that results for the parameters p describing 
the shape of the HII region are insensitive to the foreground 
model, and so have not listed additional results for models 
I, III and IV. However we find that the neutral fraction de- 
rived from observation of an HII region is sensitive to the 
foreground model assumed, and list these results in Table [2] 
Models I and II (with ni og = 2 and 3 respectively) yield 
the same neutral fraction (:ehi = 0.07 and 0.10 for evolving 
and non-evolving models respectively). However, the mod- 
els of higher order (models III and IV) remove additional 
contrast from the HII region, leading to smaller recovered 
neutral fractions. This indicates that a model which pro- 
vides for fluctuations that are of a higher order than those 
present in the true foreground will lead to removal of power 
from the inferred 21-cm signal. This in turn implies that the 
foregrounds will need to be very well understood in order 
to make quantitative measurements of the neutral fraction 
based on the observation of HII regions. 
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7 CONCLUSION 

One of the primary goals of future low-frequency telescopes 
is the detection of large quasar-generated HII regions in red- 
shifted 21-cm emission during the epoch of reionisation. Like 
all 21-cm signatures of the reionisation era, the detection of 
a quasar HII region will need to overcome the difficulties as- 
sociated with the removal of bright Galactic and extragalac- 
tic foregrounds. In this paper we have assessed the impact 
of removing a spectrally smooth foreground on the observ- 
able properties of high-redshift quasar HII regions. We have 
assumed prior removal of cxtragalactic point sources and 
considered only the unpolarised Galactic synchrotron fore- 
ground. 

The primary effect of continuum foreground removal is 
to erase contrast in the image. In particular, contributions 
to the 21-cm intensity fluctuations that have a scale length 
comparable to the frequency bandpass of the observation 
are removed by foreground subtraction. On the other hand, 
we find that this loss of contrast does not affect the ability 
of 21-cm observations to measure the size and shape of the 
HII region. We model the HII region shape as a spheroid 
described by six parameters and show that the shape re- 
covered following foreground removal agrees well with the 
shape derived directly from fitting using the input model. 

Using the recovered best-fit shape of the HII region, 
the global neutral fraction of hydrogen in the IGM could be 
measured directly from the contrast in intensity between re- 
gions that are within and beyond the HII region. However, 
we find that since foreground removal lowers the observed 
contrast between the HII region and the IGM, such a mea- 
surement of the neutral fraction would require a correction 
factor. Our results suggest that the value of this correction 
factor depends on the reionisation history. This correction 
factor would therefore need to be modeled using uncertain 
astrophysics. In addition, the measured contrast of the HII 
region, and therefore the inferred neutral fraction, is sensi- 
tive to the degree of the polynomial used for the foreground 
removal. Thus, measurement of the neutral fraction from 
quasar HII regions will require a detailed knowledge of the 
continuum foreground spectrum. 

Finally, we have considered both cases where the quasar 
is observed in an IGM which evolves slowly relative to the 
light crossing time of the HII region, and where the IGM 
evolves rapidly. The latter case is likely to be more relevant 
for observations around high-redshift quasars which are ob- 
served to be very rare at z > 6. Our results indicate that 
the evolution of the IGM will not impact negatively on the 
ability of 21-cm observations to measure the size and shape 
of quasar HII regions. 
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